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Abstract: We investigate an exactly marginal M = 1 supersymmetric deforma- 
tion of SU(A r ) M = 4 supersymmetric Yang-Mills theory discovered by Leigh and 
Strassler. We use a matrix model to compute the exact superpotential for a fur- 
ther massive deformation of the U(iV) Leigh-Strassler theory. We then show how 
the exact superpotential and eigenvalue spectrum for the SU(iV) theory follows by a 
process of integrating-in. We find that different vacua are related by an action of the 
SL(2, Z) modular group on the bare couplings of the theory extending the action of 
electric-magnetic duality away from the M = 4 theory. We perform non-trivial tests 
of the matrix model results against semiclassical field theory analysis. We also show 
that there are interesting points in parameter space where condensates can diverge 
and vacua disappear. Based on the matrix model results, we propose an exact elliptic 
superpotential to describe the theory compactified on a circle of finite radius. 



1. Introduction 



It has been known for some time that the Af = 4 supersymmetric Yang-Mills the- 
ory in four dimensions possesses exactly marginal deformations. Leigh and Strassler 
[1] identified two Af = 1 SUSY-preserving exactly marginal directions other than 
the M = 4 gauge coupling itself, giving rise to a three (complex) dimensional 
renormalization-group-fixed manifold of Af = 1 superconformal theories (SCFTs) 
of which the M = 4 fixed line is a subset. Given that strong and weak coupling on 
the Af = 4 fixed line are related by the action of SL(2, Z) on r, the gauge coupling 
of the Af = 4 theory, it is then natural to ask whether this duality extends in some 
non-trivial way over the entire A/" = 1 fixed manifold. 

In this article we answer this question for a particular Leigh-Strassler defor- 
mation of the Af = 4 theory: the so-called u q- deformation" [2] 1 . By studying a 
mass deformation of this Leigh-Strassler SCFT using a recently proposed connection 
between Af = 1 theories and matrix models [6-8], we demonstrate that Montonen- 
Olive electric-magnetic duality, extended to the full SL(2, Z) modular group, extends 
non-trivially to this space of Af = 1 superconformal theories accompanied by a well- 
defined action on the marginal parameters. Not only will we uncover the duality 
group action in these theories, but we will also provide rather powerful checks on 
the applicability of the matrix model approach towards solving these Af = 1 SUSY 
theories. (Recent checks of the matrix model proposal for other Af = 1 models have 
been performed in [3-5]). In particular we will test the superpotential and eigenvalue 
spectrum obtained from the matrix model against classical field theory results and 
find nontrivial agreement between the two. Furthermore, rather remarkably the ma- 
trix model results for the vacuum expectation values of the effective superpotential 
can be used to directly infer an exact, dynamical quantum superpotential for the 
mass deformation of the Leigh-Strassler SCFT compactified on I 3 x S 1 '. This super- 
potential written in terms of the effective fields of the 3-dimensional theory encodes 
the entire vacuum structure of the four-dimensional theory and satisfies extremely 
non-trivial checks to be discussed below. It turns out to be a natural deformation of 
the elliptic superpotential of [9] for the mass-deformed Af = 4, or Af = 1*, theory. 

Relevant deformations of SCFTs can tell us a great deal about the CFTs them- 
selves; in particular, as is well-known in the context of Af = 4 SUSY Yang-Mills and 
the so-called Af = 2 elliptic quiver models [8-11,14], the duality properties of the par- 
ent theories are inherited by the perturbed theories and manifest themselves as mod- 
ular properties of various exactly calculable observables in the perturbed theories. 
In fact, what happens in these known examples is that the holomorphic observables 

1 One such deformation has been analysed recently using matrix models in [3]. 
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in one vacuum get mapped into the corresponding observables in a different vacuum 
(realized in a different phase) under the action of SL(2, Z) with well-defined modular 
weights. The holomorphic quantities in question can often be computed exactly and 
reflect the duality symmetries of the parent theory in this simple fashion. 

The lesson to be drawn from this limited set of examples is that the consequences 
of dualities in these SCFTs are readily visible and more importantly, are readily 
calculable for certain relevant deformations of these SCFTs. With this in mind, 
rather than focusing on the SCFT itself we investigate a special mass perturbation 
of the Leigh-Strassler fixed points. This theory is the SU(iV) gauge theory with 
Af = 1 SUSY and three adjoint- valued chiral multiplets <f> + , <f>~ and $; the same 
matter content as the Af = 4 theory, but with the following classical superpotential 2 

W cl = Tr (iA<l[<f> + , &-]fi + Ml> + <f>- + fi® 2 ) , (1.1) 

where we have defined the g-commutator 

[$+ $-]f, = $+<£-e i/? / 2 - $-$+ e -W 2 (1.2) 

and where A and (3 are complex bare couplings. We also introduce the complex bare 
gauge coupling of this theory r = Amij g\ M + 6 /2n. This superpotential represents a 
deformation away from the Af = 4 point which is at A = 1, /3 — and M — \i — 0. 
It is worth noting that there are alternative relevant deformations we could have 
considered involving the operators Tr$ +2 and Tr$ -2 . These are equivalent to ( |1.1| ) 



as a deformation of the Af = 4 theory, i.e. when (3 = 0, but differ once the Leigh- 
Strassler marginal deformation is present. In particular, the resulting theories differ 
in the IR. However, as they only differ by strictly relevant operators they both flow 
to the same fixed point in the UV: namely the Leigh-Strassler SCFT. The reason for 
choosing to study the specific relevant deformation ( |1.1|) is simply that the resulting 
matrix model turns out to be exactly soluble [6, 7, 15]. 

Aside from masses for the chiral multiplets $ ± and $, our theory also involves 
two trilinear deformations, 0\ = Tr $[<l + , $~] and 2 = Tr l>{<l + , Both the 

operators 0\ and O2 are of course, marginal by power counting but only one of 
them, or more precisely, only one linear combination of these operators is an exactly 
marginal deformation of the M = 4 theory. Adding 0% to the Af = 4 Lagrangian only 
changes the coefficient in front of the Af = 4 superpotential, but it is an irrelevant 
operator even at one-loop at the Af = 4 fixed point [16]. The operator 0\ is actually 
a descendant and not a chiral primary field in the Af = 4 theory, hence its dimensions 



2 We choose a normalization in which the kinetic terms of the chiral multiplets do not have 
factors of l/S'y&r m front. In the following, hatted quantities refer to SU(iV) to distinguish them 
from U(JV). 
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are not protected. 3 The second operator O2 on the other hand is known to be exactly 
marginal at one-loop at the M = 4 point [16]. More generally, away from the M = 4 
line one should expect, for fixed Qym-, a particular linear combination of (D\ and O2 
to be exactly marginal. We have parametrized this particular linear combination via 
two complex numbers A and (3. In principle, A should be determined as a function 
of j3 and gyhi on the fixed manifold. Fortunately the specifics of this will not be 
important for us. Note that in the SU(2) theory the operator O2 vanishes identically 
and the above superpotential does not yield a marginal deformation. We will assume 
that N > 2 throughout this paper. 

The theory with masses for all the fields, as in Eq. ( |1 . 1|) , has a number of vacua, 
many of which are massive. The canonical examples are the Higgs and confining 
vacua which will play an important role in our discussion. The confining vacua 
correspond to the trivial classical solution of the F-term conditions with <f> = $ ± = 
and preserves the full SU(iV) gauge symmetry classically. At low energies the only 
classically massless fields comprise of the M = 1 gauge multiplet which confines and 
generates a mass gap. There is also a Higgs vacuum where at the classical level the 
gauge symmetry is completely broken. In addition, there are other massive vacua 
which are visible classically as solutions that leave a non-abelian gauge subgroup 
unbroken. Their classification is similar to that of the M = 1* theory [9, 10, 18]. 
What is interesting about the Leigh- Strassler deformation is that we shall find other 
vacua which are not present in the M = 1* theory itself. In tandem with this, for 
special values of the deformation parameter /3, such that e 1 ^ is a root of unity, vacua 
can disappear as a result of the condensates diverging. 

In this paper, we compute the values of the condensates (Tr$ fc ) and the effec- 
tive superpotential in all the massive vacua of the SU(iV) theory and we find that 
Montonen-Olive duality indeed relates the holomorphic condensates in different vacua 
with a special action on the couplings that we discuss below. This explicit compu- 
tation is made possible by the recent proposal of Dijkgraaf and Vafa [6] wherein the 
effective superpotential of the gauge theory is computed by the genus zero free energy 
of a holomorphic three-matrix integral [7] . The matrix model however, computes the 
superpotentials for the U(iV) gauge theory. One of our important conclusions is that 
the SU(iV) superpotential differs non-trivially from its U(iV) counterpart. (See also 
the work of [5] where a similar issue is discussed) . We show however that the former 
can be unambiguously extracted from the latter by a process of "integrating in" of 
additional fields present in the U(iV) gauge theory. This difference turns out to be 
crucial as it is the SU(iV) results that clearly exhibit Montonen-Olive duality while 



3 In the context of the AdS/CFT correspondence it is expected that these operators will get mass 
dimensions ~ (gyM^) 1 ^ 4 [^] m the strongly coupled M = 4 theory in the large- N limit. 



3 



the U(iV) results do not. For the SU(iV) theory we find in a generic (p, k) massive 
vacuum (up to inconsequential vacuum-independent additive constants) 

su ( v) = WV/xM 2 e[(p(3/2\r) 

eff 2A 2 sin / 3 0i(p/3/2|f) ' 1 ' 



with 



P TR + k iN \ \ at / n 1 i n /n 
; r iJ = r In A; p ■ q = N ; fc = 0, 1, ...,q- 1 . (1.4) 

q 7r 



The main consequence of this result is that the values of the effective superpo- 
tential (and indeed all the eigenvalues of $) in different massive vacua of the theory 
with deformation parameter (3 are mapped into one another by the action of the 
SL(2, Z) transformations on the couplings: 

ar R + b f3 2 . A 2 sin/3 

Tr — ► — ; /3 — ► —7; A sin /3 — > — , 1.5 

ctr + a ctr + a ctr + a 

ad — be = 1; a, b,c,d E Z. In particular, duality of the underlying SCFT is actually 
realized via modular transformations on a particular combination of the bare cou- 
plings t r rather than the gauge coupling r. In addition, SL(2, Z) permutes vacua of 
a theory with deformation parameter f3 provided (3 transforms with modular weight 
— 1 as above along with a specific action on A deduced from above. 

The vacuum structure and modular transformations described above can also 
be understood via an exact elliptic superpotential for the theory compactified on 
1R 3 x S 1 . The superpotential is a function of iV chiral superfields X a (a = 1, 2, . . . N, 
J2 a =i X a = 0) which parameterize the Coulomb branch of the compactified the- 
ory. These are a complex combination of the Wilson lines and dual photons of the 
compactified theory. The superpotential 

a /r 2 

W R3xS i = \ (C(*« -*b + PwiM - C(X a -X b - M/vr)) (1.6) 



a^b 



reproduces the vacua described in Eqs. (1.3) and (1.4), taking the value (1.3) de- 
scribed above. Here, ((z) is the Weierstrass-C function for the torus with half-periods 
i>j\ = m and 0U2 = ^tr, i.e. with complex structure Tr. This superpotential also 
predicts new SL(2, Z)-invariant vacua, not present in the Af = 1* theory, whose 
existence is confirmed by the classical analysis of the four- dimensional field theory. 

We also remark that the properties of (3 under modular transformations are 
similar to those of individual gauge couplings in elliptic quiver theories [14,19]. In 
addition, at non-zero (3 the matrix model solution naturally involves a torus with 
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two marked points with (3 parameterizing the separation. These suggest a deeper 
connection between the Leigh- Strassler and quiver SCFTs. 

The layout of this paper is as follows: in Section 2 we show how one can relate the 
superpotentials for the U(iV) and SU(iV) theories by a process of integrating-in the 
trace part of the fields. This is very important because the matrix model approach 
yields the superpotential in the U(iV) theory, whereas, we are primarily interested 
in the SU(iV) theory. In Section 3, following Dijkgraaf and Vafa, we show how a 
matrix model can be used to find the superpotential in the confining vacua of the 
U(iV) theory and therefore by implication in SU(iV) theory. In particular, we find the 
spectrum of eigenvalues of the adjoint- valued field $ from which all the holomorphic 
condensates can be calculated. Section 4 is devoted to an analysis of the vacuum 
structure of the mass deformed Leigh-Strassler theory from the point-of-view of the 
tree-level superpotential. In particular, this leads to expressions for the classical 
limit of the condensates in each of the massive vacua. In addition, for the case of 
gauge group SU(3) we find new vacua that are not present in the Af = 1* theory. 
We also show how some of the massive vacua can disappear when the coupling (3 
takes particular values. In Section 5, we return to the matrix model and consider 
multi-cut solutions that describe all the massive vacua. We then go on to how the 
structure of vacua can be used to deduce the action of the SL(2, Z) duality group on 
the Leigh-Strassler theory itself. The final Section briefly reports on how the results 
from the matrix model can be used to deduce the exact elliptic superpotential (1.6) 
for the theory compactified on a circle to three dimensions, generalizing the one for 
the Af = 1* theory constructed in [9]. 



2. Relation between the ~U(N) and SU(iV) theories 



As we have explained one of our aims is to compute the quantum superpotential and 
a certain set of holomorphic observables in the massive vacua of the mass deformed 
SU(iV) Leigh-Strassler SCFT with Af = 1 SUSY. For this purpose, we would like to 
employ matrix model techniques that have been proposed recently by Dijkgraaf and 
Vafa [6]. However, already there is a subtlety. A direct application of the matrix 
model approach will solve the U(iV) theory since the DV proposal relates the effective 
superpotentials for U(iV) Af = 1 gauge theories to the planar diagram expansion of 
corresponding matrix models. So the first question that we must address is how are 
the SU(iV) results related to those of the U(iV) gauge theory? We will now show that 
there is a very specific relation between the effective superpotential of the SU(iV) 
theory with classical superpotential Eq. (|1.1|), and the effective superpotential of the 
U(iV) theory with the same classical superpotential. This relation will eventually 
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allow us to extract the SU(iV) results from the matrix model for the U(iV) gauge 
theory. 

2.1 The U(iV) gauge theory 

The Leigh-Strassler theory with SU(iV) gauge group differs non-trivially from its 
U(iV) counterpart. The U(iV) theory contains additional neutral chiral multiplets 
that couple to the chiral superfields transforming in the adjoint representation of 
the SU(iV) C U(N). These interactions can modify the superpotential and other 
holomorphic observables of the SU(iV) C U(iV) theory. The way this happens is 
clearly seen from the point-of-view of the U(iV) theory. Let us begin by considering 
the U(iV) version of the theory, namely, the M = 1 SUSY gauge theory with tree 
level superpotential 

W V(N) = Tr ^ A$ ^+ ) $ -]^ + M $+<J>- + ^$2) ^ 

where $ ± and $ are the fields in the adjoint of U(iV). Now the fields $ ± and $ can 
be naturally split into their traceless and trace parts: 

$ = <f> + a; $ ± = <l ± + a ± ; Tr<f> = 0; Tr^^O. (2.2) 

Here, $ ± and $ are traceless and so transform in the adjoint representation of 
SU(iV) C XJ(N) while a ± = Tr^/N and a = Tr $/iV are neutral. Rewritten in 
terms of these variables the tree level superpotential fl2.1|) for the U(iV) theory is 

W V(N) = N ( M _ 2Aasinf)a+a- + Nfia 2 + W* v (JV) (a, a ± ) (2.3) 

where 

H/J U(JV) (a,a ± ) =Tr fi\^ + , $~]^ + [M - 2Aa sin |) $ + <l _ + yU<l 2 

x (2-4) 
-2Asinf a"$$ + -2Asinf a + $$-j . 

The main point here is that the neutral trace fields a and a have the effect 
of modifying the couplings of the SU(iV) fields. For example, the mass M has been 
renormalized to M— 2Xa sin |. In addition, there are new bilinears in ^ and $ whose 
couplings 4 depend on a, a 1 * 1 . (Notice that when (3 = 0, that is for the Af = 1* theory, 
these modifications disappear and there is no real difference between SU(iV) and 

4 Strictly speaking, of course, these "couplings" are actually chiral superfields. But from the 
point of view the SXJ(N) sub-sector of the U(iV) theory, these neutral chiral superfields do appear 
like couplings that have been elevated to chiral superfields otherwise known as "spurions" . 
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U(iV).) In fact, we show below that because of the symmetries of the theory the VEVs 
of a + and a~ are forced to be zero self-consistently in the full theory. Thus the only 
effect of the neutral fields is to modify the mass term M — > M' = M — 2a\ sin ~. One 
might then suspect that answers for pure SU(iV) theory of Eq. (|1.1|) may be obtained 
by simply rescaling the U(iV) results by appropriate powers of M' = M — 2a\ sin |. 
We will now see that this is almost correct. 

First of all, the fields a, a ± are blind to all gauge interactions while $, $ ± ex- 
perience only SU(iV) gauge interactions, the remaining U(l) gauge multiplet being 
decoupled from everything else. Hence in a vacuum where the SU(iV) fields are 
rendered massive and the gauge interactions generate an effective superpotential for 
that sector of the theory, we may readily write the effective superpotential for the 
U(iV) gauge theory as 

W^ N) =N(M- 2Xa sin § )a + a~ + Nfxa 2 + W^ (N) (a, a ± ) . (2.5) 

Note that the chiral fields a, a 1 * 1 appear as parameters or couplings for the SU(iV) 
sub-sector, however they are actually dynamical variables in the full theory and 
their values must be determined by extremizing the full effective superpotential. 
This procedure can actually be implemented formally by first noting that the SU(iV) 
sub-sector, from Eq. ( |2.4| ), has certain abelian symmetries. Let us define for the 
sake of convenience M' = M — 2Aasin|. It is then sufficient to consider the fol- 
lowing discrete symmetries: (i) (<l )± ,a ± ) — > — ($ =t , a^); (ii) ($ + , $, a + , M', fx) — > 
(-i<& + ,i<&,-ia + ,iM',-fx); and (iii) (l>~, <f>, a~, M', fx) — > (-i < l _ , i<&, —ia~, iM', -fx). 
The only possible form for the SU(iV) effective superpotential, consistent with these 
symmetries, analytic in the parameters and which has mass dimension three, is 

W^ (N) = fxM' 2 F(a + aT IfiM') , (2.6) 

for some unknown function F. Plugging this back into Eq. ( ^.5[ ) for the effective U(iV) 
superpotential, and imposing the F-term conditions by extremizing with respect 
a, a^, the only solution that generates a nontrivial effective superpotential is the one 
where 

(a±> = 0; and ^ J (a) = wT'- (2.7) 
2A sin | 

This tells us that in a vacuum of the U(iV) gauge theory, the trace of the adjoint 
scalar $ must be related to the value of the effective superpotential of the U(iV) 
theory in that vacuum, precisely according to the above equation. This already 
constitutes a non-trivial prediction for the Dijkgraaf-Vafa approach for solving the 
U(iV) gauge theory, one that our results from the matrix model must satisfy. 
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2.2 The relation between the U(N) and SU(iV) superpotentials 



The symmetry arguments above also imply that the effective superpotential of the 
SU(iV) theory with a classical action as in Eq. fll.lQ must be 

W^ (N) = I^M 2 F(0) . (2.8) 

But we can easily determine the unknown function F in terms of the U(iV) super- 



potential using Eqs. (|2.5|) , (|2.6|) and ( |2.7|) and we find: 



W SV{N) = (2 9) 

eff , 4A 2 sin 2 /3/2 u/ U(Af) ' 

1 NfiM 2 W eS 

Anticipating the matrix model results of the following section, we will write this in 
a slightly different way that may turn out to be more illuminating: 



ci 

4A 2 sin 2 g/2 T TrU(iV) 



W SU(N) _ M 



cff 



M' 2 ((a)) 



eff 4A 2 sin 2 | 



NuM 2 

■ (2.10) 

4A 2 sin 2 f 1 ; 



In this form, it is apparent that the effective superpotentials of the U(iV) and SU(iV) 
gauge theories are indeed related by the replacement M'((a)) — > M after subtracting 
off certain additive constants. This form of the relation will turn out to be quite 
suggestive and useful when we discuss the matrix model results. Importantly, this 
simple relation tells us how to extract the SU(iV) answer from the U(iV) result which 
the matrix model naturally computes. 



3. The U(iV) theory from the matrix model 

According to the Dijkgraaf-Vafa proposal [6], in any given vacuum, the effective su- 
perpotential of the U(iV) Af = 1 gauge theory with classical superpotential Eq. ( |1.1| ) 
is computed in terms of the planar diagram expansion of the three-matrix model 
partition function expanded around that vacuum 

Z = J [d$ + ] [d$~] [d<$>] exp -^Tr (zA$[$ + , $~]^ + M$ + $~ + /i$ 2 ) . (3.1) 

We use the same notation for the matrix fields and the associated superfields in the 
U(iV) theory. In the matrix model, unlike the field theory, one takes $ + = ($ - )^ 
with the fluctuations of $ around the saddle point to be Hermitian. This matrix 
model has been actually solved in a different context [15] in the large- N limit. First 
one integrates out <I >± and performs a field rescaling $ — ► $/ A to get 

Z = X- N2 f\d*] exp-^VTr$ 2 /A 2 

J 1 1 |det(Ml®l-ze-^/ 2 $®l + ^/ 2 l®<l>)| ' 1 ' 



S 



Now we can follow [15] to obtain the saddle-point equation in the large- N limit. For 
completeness we will now follow the steps required. 



Let {4>i} denote the eigenvalues of the matrix $. Changing the integration 
variables in the matrix integral by going to the eigenvalue basis introduces a Jacobian: 
the famous Van der Monde determinant which leads to a repulsive force between 
the eigenvalues. The second step is a variable change that will eventually yield a 
simplified form for the saddle-point equation: 

M 

<f>i = 'Me 5 * + — —g . (3.3) 
2 sin I 

In terms of these variables the classical potential for the matrix model eigenvalues 
takes the following form 

^E^ + J^, (3.4, 
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where we have defined 

In addition to this classical potential, the eigenvalues 5i also experience pairwise 
effective interactions induced both by the Van der Monde determinant and the de- 
terminant resulting from integrating out $ ± in Eq. (|3.2j ). The eigenvalues 5i are 
naturally defined on the complex- z plane with the identification z ~ z + 2ni, i.e. a 
cylinder. 



3.1 Solving the large-iV matrix model for the confining vacuum 

In the large- N limit, the eigenvalues form a continuum and condense onto cuts in 
the complex plane. On can think of these cuts as arising from a quantum smearing- 
out of the classical eigenvalues of $. For the confining vacuum all the classical 
eigenvalues are degenerate, since $ cx 1, and so we expect a solution in the matrix 
model involving a single cut. Multi-cut solutions will be discussed later. The extent 
of the cut and the matrix model density of eigenvalues p(5) can be determined from 
the saddle-point equation in terms of the parameters of the classical potential and 
the 't Hooft coupling of the matrix model S = g s N . The saddle-point equation is 
most conveniently formulated after defining the resolvent function: 

u>{z) = \ jT 'd* jffL , Se[a,b], J'p(5)d5=l. (3.6) 
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Figure 1: The region over which the function G(z) is defined with its two cuts. The lines 
Imz = ±iir are identified. 

This function is analytic in z and its only singularity is along a branch cut extending 
between [a, b]. The matrix model spectral density p(S) is equal to the discontinuity 
across the cut 



u(8 + ie) — u(5 — ie) = —2iiip(5) 



5 e [a, b] 



(3.7) 



The saddle-point equation expresses the zero force condition on a test eigenvalue in 
the presence of the large- N distribution of eigenvalues along the cut. For the matrix 
model that we are studying, this condition is best written in terms of the function 
G(z) defined as 

G(z) = U(z) + iS(uj{z + %\) - oo(z - if)) (3.8) 
and U(z) is a polynomial in e z such that 



V'( z ) = -i(U(z + il)-U(z-il)) 



(3.9) 



From this definition, one may easily deduce that G(z) is an analytic function on 
the cylinder and that in the interval |Imz| < it it has two cuts [a + if, b + if ] and 
[a — if, b — if ]. This is illustrated in Figure 1. In terms of G(z), the matrix model 
saddle-point equation is 



G(5 + zf ±ie) = G{8 



5 e [a, b] 



(3.10) 



The saddle-point equations actually provide a gluing condition on the cylinder with 
two branch cuts. In particular, this defines a torus with two marked points. G(z) is 



then uniquely specified by gluing condition Eq. ( 3.10 ) and asymptotic behaviour at 
large z 



lim G{z) 



/iM 2 



1 



sin (3 



:e 2z - 



1 



2 sin 2 f 



0(e 



(3.11) 
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Figure 2: Uniformizing map from the z-plane to the torus Ef. 

which is a consequence of Eqs. (|3.8|) and ( |3.9| ). With this data we can determine 
G(z) in a suitable parametrization so that we can then implement the DV proposal 
and obtain the U(iV) gauge theory superpotential. 



3.2 The r-torus and the elliptic parametrization 



As we stated earlier, the two-cut cylinder along with the asymptotic properties of 
G(z) in Eq. ( pMlp uniquely specifies a torus E f with complex structure parameter 
f. As shown in Figure 2, the contour Ca enclosing one of the cuts [a + i|, b + i^} 
anticlockwise maps to the A-cycle of the torus while the contour Cb joining the two 
cuts maps to the S-cycle. We will think of this torus as a quotient of the complex 
w-plane (u being an auxiliary variable) by a lattice with periods 2o)i and 2o>2 and 
with complex structure f = lj 2 /^i- 

^From the field theory view-point, the appearance of this torus is not entirely 
obvious except for the fact that the gauge theory we are studying is connected to the 
Af = 4 theory. In the Af = 1* case discussed in [7,8], one might have expected the 
torus to appear due to the relation with the Seiberg-Witten curve of the associated 
Af = 2* theory. For the Leigh- Strassler SCFT the fact that the torus structure is 
also present is one of the most important conclusions of this paper: as we shall see 
it leads to an understanding of how SL(2,Z) acts on the underlying SCFT. 

The map z(u) from the w-plane to the z-plane is specified by the requirements 
that going around Ca (the A-cycle of the torus) returns z to its original value, while 
traversing the contour Cb (the S-cycle of the torus) causes z to jump by an amount 
i/3, which is the distance between the two cuts in the z-plane. Both these operations 
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leave G unchanged which is therefore an elliptic function on the u-plane. Thus 



A - cycle : z[u + 2u>i) = z(u) ; G(z(u + 2&i)) = G(z(u)) , (3.12a) 

B -cycle: z(u + 2cD 2 ) = z(u) + i(3 ; G(z(u + 2Cj 2 )) = G(z(u)) . (3.12b) 

This determines the following unique map z(u) from the w-plane to the two-cut 
cylinder: 

exp z(u) - -J- . mm^u^-p/Mf) 

The fact that G(2;(m)) is an elliptic function on the w-plane along with its asymptotic 
behaviour at large z ( |3.11|) (corresponding to a pole in the w-plane) leads uniquely 
to 

GWm)) = iSSf ' twS w " + ^' 8/2 ^ ~ ' (3 ' 14) 

Here, p(u) is the Weierstrass function, an elliptic function defined on Ef, and 9i(u\t) 
is the Jacobi theta function of the first kind (see [20] and Appendix A for more 
details). As expected, G(z(u)) is an elliptic function of u with a second order pole in 
the li-plane; z(u) on the other hand is only quasi-elliptic. For a different choice of the 
deformation potential the asymptotic behaviour of G(z) will change and G(z(u)) will 
end up being an elliptic function of a different order, but significantly z(u) will remain 
unchanged. There are two special points on the torus Ef at u — ±/3d)i/27r which map 
to the points z = =Foo. Having determined G(z) in the elliptic parametrization we 
can now implement the DV proposal to compute the superpotential in the confining 
vacuum of the U(iV) gauge theory. 



3.3 Superpotential of the U(N) theory via Dijkgraaf-Vafa 

According the Dijkgraaf-Vafa proposal, the gluino condensate of the gauge theory 
gets identified with the 't Hooft coupling S of the matrix model. From Eqs. ( |3.6| ) 
and (|3.8|) , the integral of G(z) around one of its cuts, say Ca, is equal to —2^5*, so we 
can directly compute the gluino condensate of the U(iV) gauge theory in the elliptic 
parametrization: 

U A = 2mS = -i I G(z(u))^P-du . (3.15) 
J A du 

The second ingredient required to determine the QFT superpotential is the variation 
in the genus zero free energy of the matrix model in transporting a test eigenvalue 
from infinity to the endpoint of the cut. This is obtained by integrating the force 
on a test eigenvalue from infinity to the endpoint of the cut and can be expressed 
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as II B = — i j B G(z)dz. There is a slight subtlety here as there is a "zero-point" en- 
ergy associated with additive constants in the matrix model potential. First note that 
there is a multiplicative factor of A"^ 2 in the partition function ( |3.2|) and an additive 
constant in the classical scalar potential ( |3.4| ). Both these appear as S'-dependent con- 
tributions to the matrix model genus zero free energy defined as F /g 2 s = N 2 F /S 2 , 
so that 

^ = U B -2S\n\+ ; U B = -i[ G(z(u))^-du. (3.16) 

oS 4A 2 sm 2 f J B du 

The effective superpotential of the mass deformation of the U(iV) Leigh- Strassler 
theory in its confining vacuum is obtained by extremizing the following expression 
with respect to S: 

WT> = * §| - 2mrS; ±W™ = ° ■ (*■") 

where r is the bare coupling of the theory. The ^-independent constant in Eq. ( |3.16i ) 
appears as a vacuum-independent additive constant in the superpotential. On the 
other hand, the term linear in S in Eq. ( |3.16|) actually has the effect of renormalizing 
the bare coupling r so that 

iN 

r — > t r = t In A . (3.18) 

7T 



Both Ua and IT^ are easy to evaluate in our elliptic parametrization: 

n, = ; U B = r ^ - h(f) , (3.19) 

dr dr 



where 



_ ^cosf e^/m 

() ~4A 2 sin 3 f W/m' [ } 



It follows straightforwardly that 

OS eff N 



4^T ) = => r= T 4 (3.21) 



so that 

WT = -NMr R /N) + . (3.22) 

Finally we have, therefore, the quantum superpotential of the U(iV) gauge theory in 
its confining vacuum, 

U(JV) iV>M 2 cos f 6 1 (/3/2\r R /N) iV/iM 2 

eff 4A 2 sin 3 | ' 6[((3/2\t r /N) + 4A 2 sin 2 f ' 1 ' 
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with tr defined in ( |3.18|) . This result satisfies the first obvious check: in the (3 — > 
limit with A = 1 it reproduces the superpotential of the N = 1* theory. The 
superpotential in the other N—l confining vacua are obtained by repeated T-dualities 
tr — > tr+1. Note that the T-dualities in question naturally act on tr rather than the 
bare coupling r. Similarly the result for the superpotential has a modular parameter 
tr/N and so we expect that any modular properties we uncover should naturally be 
a consequence of the modular group acting on tr the "renormalized" coupling. 



3.4 Eigenvalues of $ in the U(N) theory and a test 

One of the main results of our earlier work [7] in the context of M = 1* theory 
was that the map z(u) from the w-plane to the z-plane, when evaluated along the 
A-cycle of the f-torus (which maps to the upper cut on the z-plane), precisely yields 
the eigenvalues of the adjoint scalar $ of the gauge theory in its confining vacuum, 
distributed uniformly along this cycle. 5 We were able to prove this assertion for a 
class of deformations of the M = 4 theory involving masses for two of the adjoint 
chiral superfields and arbitrary polynomials in the third adjoint field $. We will 
not pursue this generalized analysis in the present case, but we will present strong 
evidence that the same relation with the field theory eigenvalues holds under the 
Leigh-Strassler deformation. 

For the Leigh-Strassler theory, we first recall that the 5i are valued in the z-plane 
and these are in turn related to the eigenvalues of the $ matrix through Eq. ( |3.3| ). 
Taking into account the fact that the function G(z) is defined in Eq. (|3.8|) in terms 



of the resolvent functions u(z ± i|) where the positions of the cuts of the latter 
have been displaced by ±i| , we conjecture that the function — M exp (z(u) — if)/A 
evaluated along the A-cycle gives the eigenvalues of $ in the U(iV) theory: 



<f)(x) 



2A sin f 



x e [- 



.1 li 

2' 21 



p VMtrJN) 6 3 (ttx-[3/4\tr/N) 

fOS — -r- 

2 6> x {P/2\tr/N) e 3 (nx + {3/4\T R /N) + 

(3.24) 

Note that the factor A takes into account the field rescaling prior to Eq. ( |3.3j ) and we 
also incorporate the additive shift of Mj (2 A sin |) involved in the variable change in 
Eq. (|3.3j). This is of course a continuous distribution of eigenvalues appropriate for 
the U(iV) theory in the large- iV limit. But the finite iV answer is obtained by simply 
replacing the continuous index x by a discrete index. For this to be legitimate, it 
is essential that the field theory eigenvalues be uniformly distributed in the interval 



5 In contrast the matrix model eigenvalues are not distributed uniformly along this cycle signalling 
the crucial difference between the spectral densities of the two systems. Nevertheless, the matrix 
model provides a means to extract the field theory distributions in a simple way. 
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x G [— |, |] as we have indeed pointed out earlier. The difference between finite- iV 
and large- N manifests itself in gauge-invariant observables as vacuum-independent 
operator mixing ambiguities that vanish in the large- N limit [7]. 



Our results for the eigenvalues in the confining vacuum of the U(iV) theory 
satisfy several checks. Firstly in the classical limit r R — > zoo, using the properties of 
the Jacobi theta functions (see Appendix A), Eq. (|3.24j) vanishes as expected, since 
the confining vacuum should correspond to the trivial classical solution. Secondly, 
in the (3 — > limit with A = 1 the eigenvalues of M = 1* theory in its confining 
vacuum reproduce the expressions in [7,21]. Two further nontrivial tests involve the 
computation of (Tr$) and (Tr$ 2 ) in the U(iV) theory using the eigenvalues above. 
With a large- N distribution, the trace reduces to an integral and we find 



<Tr$> = N f' 2 dx<P(x) = ^^iwT ) , (3.25a) 

J —1/2 V M 

(Tr $ 2 ) = N f' 2 dx <j) 2 (x) = - W^ {N) , (3.25b) 

7-1/2 V 

where W e g is the expression obtained in Eq. (3.23). We show how to compute the 



associated integrals explicitly in Appendix B. The first (|3.25a|) demonstrates that the 
trace of $ calculated using the matrix model results satisfies precisely the relation 
Eq. (|2.7|) predicted on general grounds from field theory. The second integral is also 
exactly in accord with field theory expectations since from our general arguments in 
Section 2, we have oc \i and therefore differentiating with respect to \i gives 

(Tr$ 2 ) = W^ s {N) /jj precisely (p5TJ). 



In summary, the eigenvalues of the adjoint scalar $ in the confining vacuum of 
the mass deformed U(iV) Leigh-Strassler theory at large iV are given by Eq. ( |3.24j ). 
The eigenvalues of the U(iV) theory at finite iV are obtained by substituting the 
continuous index x with an appropriate discrete index. 



3.5 The SU(iV) superpotential 



We argued in Section 2 that the superpotentials of the U(iV) and SU(iV) theories 
differ non-trivially, but are nevertheless simply related by Eq. (|2.10|) . The main point 
there was that the U(iV) theory allows a non-zero (Tr $) = N(a) ^ to develop in 
the confining vacuum of the quantum theory. This leads to what is effectively a 
renormalization of the mass parameter M — > M' = M — 2(a)Asin|. According 
to the result of our field theory analysis of Section 2, the SU(iV) superpotential 
may be recovered by a rescaling of the U(iV) result by M 2 /M' 2 , after subtracting 
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off a vacuum independent constant. It is quite interesting to note that precisely the 
same vacuum independent constant makes an appearance in the U(iV) superpotential 
( |3.23|) computed in the Dijkgraaf-Vafa approach. 

Using our results for the U(iV) superpotential from the previous section and 
Eq. ( p.lO| ), we thus find that in the confining vacuum the SU(iV) theory develops the 
superpotential, 

su W N^M 2 9[(P/2\t r /N) _ jV/xM 2 

eff ~~2A 2 sin/? 9 1 (P/2\r R /N) 4A 2 sinf' l ' J 

where the renormalized coupling is defined in ( |3.18| ). This expression for the SU(iV) 
superpotential has the correct /3 — limit in that it reproduces the Af = 1* super- 
potential exactly. 



3.6 SU(iV) eigenvalues and condensates 

The eigenvalues <pi of the adjoint scalar <f> in the confining vacuum of the SU(iV) 
theory are related to the U(iV) eigenvalues by a simple extension of our arguments 
in Section 2 as 



4>(x) 



M 
~M' 



M 



[x - 



2A sin f 



+ 



M 



2A sin f 



(3.27) 



This can be understood simply as a rescaling of the U(iV) spectrum of eigenvalues 
by M/M'((a)) after subtracting out a constant and is in line with the arguments of 
Section 2 for extracting the SU(iV) superpotential. 



Evaluating this explicitly yields 

6'Mt r /N) 6 3 (irx - /3/4\t r /N) 



4>{x) 



M 
2A 



e 1 {(5/2\T R /N) 6 3 (nx + P/4\r R /N) 



+ 



sin 



x e - 



1 l \ 

2' 2i 



(3.28) 

and (Tr$ 2 ) = 

W^'/fi. Once again we remark that with a continuous label x these are really the 
eigenvalues of the large- N gauge theory but finite- N results are obtained by simply 
replacing x with a discrete index. 



We show in Appendix B that with these eigenvalues (Tr<J>) 

rSU(iV) 



Now that we have the complete distribution of eigenvalues of $ in the confining 
vacuum of the SU(iV) theory, every condensate of the form (Tr$ fc ) is automatically 
computed via the integrals 



(Tr = N 



1/2 



dx 4> k {x) 



(3.29) 



-1/2 
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We will not undertake the exercise of computing these but, if desired, they can be 
computed following the methods in Appendix B. 

3.7 Summary 

In this section we have computed the value of the effective superpotential and the 
spectrum of eigenvalues of $ in the confining vacuum of the SU(iV) theory. We 
did this by obtaining the U(iV) superpotential via a matrix model and subsequently 
extracting the SU(iV) results by a process of integrating-in 6 discussed in Section 2. 
The results of the computation have satisfied several general field theory checks. In 
addition, the emergence of a quasi-modular function (9[/9i) as the superpotential 
of the SU(iV) theory is strongly suggestive of Montonen-Olive duality being realized 
in this theory. (In contrast, the U(iV) superpotential being basically the reciprocal 
of the SU(iV) result, transforms in a complicated way under the modular group.) 
In order to establish electric-magnetic duality we need to be able to compute the 
corresponding quantities in other massive vacua of the theory (most notably the 
Higgs vacuum which should be S-dual to the confining phase). We will address this 
in Sections 4 and 5 which will be devoted to understanding the classical and quantum 
vacuum structure respectively. 



4. Vacuum structure of deformed SJJ(N) Leigh- Strassler the- 
ory 

In this section we exhibit the classical vacuum structure of the SU(iV) theory. Im- 
portantly these classical configurations will serve as the starting point for analysing 
other saddle points of the matrix model and thence obtaining the effective quantum 
superpotential in all the massive vacua. Interestingly the classical solutions (the 
eigenvalues) for the Higgs vacuum that we find here will be shown to be related by 
S'-duality to the quantum eigenvalue spectrum of the confining vacuum that we have 
already obtained in the previous section. 

The vacuum structure of the mass deformed SU(iV) Leigh-Strassler theory turns 
out to be subtly different from that of the N = 1* theory: there are new vacua and 

6 We remark that one might have expected to obtain the SU(iV) superpotential and eigenvalues 
directly from the matrix model by imposing a constraint on the matrices via a Lagrange multiplier. 
However this turns out not to work in any simple way, the primary reason being, as pointed out in [7] , 
that the matrix model eigenvalue spectrum is completely different from the field theory spectrum 
and tracelessness of the latter translates into a rather complicated constraint on the former. 
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old vacua can disappear. In order to describe the vacuum structure, we start with 
the tree-level superpotential (|1.1|). At the classical level vacua are determined by 
solving the F- and -D-flatness conditions modulo gauge transformations. As usual, 
if one is interested in holomorphic information it is sufficient to solve the F-flatness 
conditions modulo complex gauge transformations. These conditions follow from 
extremizing the tree- level superpotential ( |1 . 1| ) : ' 

P[$ + ,$-bj = -i^<£ , P [$, $+] = , P [<£>",$] = -i— $~ . (4.1) 

A A A 

In the above, P is the projector onto the traceless part which must be imposed 
since we are dealing with the SU(iV), rather than the U(N), theory. Complex gauge 
transformation can be used to diagonalize $. 

In the limit (3 = 0, the theory becomes the basic M = 1* theory and the condi- 
tions ( [Lip reduce to the Lie algebra of SU(2). Let us briefly recall the solutions in 
that case. Each solution is associated to a representation of SU(2) of dimension N 
where $ and ^ are, up to a trivial re-scaling, J3 and J±, respectively. Each repre- 
sentation of SU(2) of dimension N is associated to a partition iV — > rii + ■ • • + n q , 
where each integer rij is the dimension of each irreducible component. In general the 
unbroken gauge group has abelian factors and the phase is Coulomb, or massless. 
The massive vacua are associated to the the equi-partitions N = q ■ p for which the 
unbroken gauge group is SU(g). In such a vacuum, the low energy description is pure 
Af = 1 with gauge group SU(g) which has q inequivalent quantum vacua. Hence, 
there are J2 q \N 1 vacua with a mass gap (the sum being over the positive integer 
divisors of N) . 

Now we consider the case when j3 7^ 0. Each of the vacua of the M = 1* theory 
remains a vacuum of the deformed theory, at least for generic values of f3. These 
deformed solutions have the following form. Firstly, the block form of the solutions 
associated to the partition N — > n% + • • • + n q is preserved. Let us concentrate on 
one of the blocks of dimension n. For /3 — 0, $ and ^ are equal to J3 and J±, up 
to trivial re-scaling, of the irreducible SU(2) representation of dimension n: 

P = * = *±=^ ± . (4.2) 

When j3 7^ 0, the non-zero elements of $ ± have the same pattern as J± (so only 
the elements just above and below the diagonal, respectively are non-zero) but the 
numerical values are changed. The elements of diagonal matrix $ in the block in 



7 Note that the F-term equations are rather similar to the "g-deformed" SU(2) algebra. They 
are similar to, but not the same as the equations encountered in [2]. The authors of [2] considered 
a slightly different mass deformation and a U(iV) gauge group. 
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question are found to be 



M ( 1 n 



- 



3 2AV sin f sinf 



-i/3(n+l-2j)/2 



j = l,...,n. (4.3) 



Notice that when (3 — > 0, in this block, becomes equal to — (M/A)J 3 as expected. 

For generic (3, these deformed vacua are massive/massless according to their type 
when (3 = 0. However, for specific values of (3 vacua can migrate off to infinity in field 
space and hence cease to be genuine vacua of the field theory. This occurs when the 
eigenvalues of $ diverge. It is easy to see that this will occur in a block of dimension 
n when the deformation parameter 

9-jrk 

13 = — , keZ. (4.4) 
n 

In one of the N/p massive vacua with unbroken SU(N/p) symmetry the eigen- 
values of $ are given by (|4.3j ) with n = p and where each eigenvalue has degeneracy 
N/p. From this spectrum, we can calculate the classical values for the condensates; 
for instance, for the quadratic condensate: 

, -o, M 2 N ( tanf \ , . 

Tr $ 2 = jt (p \ - 1 ) . 4.5 

V 1 4A 2 sin 2 f Vtan f / V ; 

Notice that this condensate diverges when (3 = 2irk/p, k G Z in accord with 



4.1 New classical branches 



As well as the original massive vacua, there are also new vacua which are not present 
when (3 = 0. We illustrate this is in the simplest case where the gauge group is 
SU(3). In this case there is a new vacuum, more precisely a moduli space of vacua, 
in which $ and $ ± are all, up to complex gauge transformations, diagonal: 

2 A sin 

a^r 1 ) , (4.6) 

Here, p and £ are two complex moduli of the solution and Xi, i = 1,2,3 are the three 
distinct roots of the cubic 

x 3 -3x + 2-p£ = . (4.7) 



V>M/2 



2 

A sin I 



pdiag((l-x 1 )- 1 ,(l-x 2 )- 1 ,(l 



A sin | 

9 

£diag((l-xi 



-i 



;i-* 2 )-\(i- 
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Note that the tracelessness of $ and $ ± is guaranteed by the form of the cubic (|4.7|) . 
Clearly these solutions diverge to infinity when (3 — > 0. However, the moduli space 
of the solutions also has additional branches which occur whenever 

e i ? /2 x l -e- i P /2 x J =2tsm§ . (4.8) 

In this case, an off-diagonal component of either of ($ ± )ij can become non-zero. The 
quadratic condensate in this space of vacua is 

3 /If 2 

SU(3): (Tr$ 2 )= . (4.9) 

2A 2 sin 2 f 

Clearly this diverges when (3 = 0. We will see that these new vacua are important 
for checking a conjecture that we make later. 



5. Electric-magnetic duality in the Leigh- Strassler Deforma- 
tion 



The question we now address is what becomes of the 5-duality of the Af = 4 theory 
under the Leigh- Strassler deformation? Specifically we are considering a deformation 
of the A/" = 4 Lagrangian of the form; 

C = C N=A + J d 2 9AW (5.1) 

where the deforming superpotential can be written as 

AW = i(Acosf - 1)0! - Asinf£> 2 + MTrl> + <f>- + fiTr® 2 (5.2) 

where, as in Section 1, O x = Tr <f>[<l + , $~] and 2 = Tr ${<l + , $~}. At the A/" = 4 
point, the operators appearing in the deforming superpotential transform under S- 
duality in a definite way. Adding these operators to the superpotential therefore 
breaks modular invariance. However, as is standard when dealing with broken sym- 
metries, we should be able to restore invariance by allowing the corresponding cou- 
plings to transform in an appropriate way. This strategy was implemented in detail 
for the Af = 1* case /3 = in [13]. We will now review this briefly and extend this 
discussion to the marginal deformation. 

Using the conventions described in [13] , 8 the superspace measure d 2 9 has holo- 
morphic modular weight —2. Hence for modular invariance of the Lagrangian, the 

8 See the final paragraph of Section 3 of this reference. 
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superpotential must have weight +2. In the M = 4 theory, the mass operators 
Tr $ + $~ and Tr $ 2 are both part of the multiplet of chiral primary operators which 
transform as second rank symmetric traceless tensors under the SO (6) R-symmetry. 
Their modular weights can be deduced via the AdS/CFT correspondence by com- 
parison with the dual fields in IIB supergravity [12]. In fact they have holomorphic 
modular weight +2 and the Lagrangian will be invariant if we assign the mass pa- 
rameters M and \i weight zero. Although there is no explicit breaking of modular 
invariance, it is important to remember that the theory has various vacuum states 
in which these mass operators acquire expectation values, thereby 'spontaneously' 
breaking 5-duality. In other words, modular transformations map the observables 
in one vacuum to those in another. Indeed, this property is manifest in the exact 
results of [11, 13]. 

Similar arguments should apply to the Leigh-Strassler theory at least near the 
M = 4 point where the modular properties of the deforming operators are known. 
If we set A = 1 and work to linear order in (3, the deformation only contains the 
cubic operator Oi- This operator is also part of an M = 4 chiral primary multiplet 
consisting, in this case, of third rank symmetric traceless tensors of 50(6). In the 
conventions of [13] it has modular weight +3. To restore invariance under 5-duality 
we should therefore assign the parameter (3 holomorphic modular weight —1. In Sec- 
tion 5.2, we will find the appropriate generalization of this transformation property 
throughout the parameter space, by using the matrix model to find the exact super- 
potential in each of the massive vacua. In Section 5.1 we see what can be deduced 
by looking at the Higgs and confining vacua alone. 



5.1 A first look at 5-duality 



Let us concentrate on the Higgs and confining vacua. In the M = 1* limit, it is 
well known that each observable in the Higgs vacuum is simply related by 5-duality, 
t — ► — 1/t, to the corresponding observable in the confining vacuum. This is an 
explicit realization of the Montonen- Olive electric-magnetic duality of the M = 4 
theory. The question is whether a similar duality persists under the Leigh-Strassler 
deformation. 



If 5-duality is realized in this theory it should exchange the Higgs and confining 
vacua. Since we already have the eigenvalues in the confining vacuum of the SU(iV) 
theory a simple and powerful test emerges: use 5-duality to deduce the eigenvalues 
in the Higgs vacuum and then take the classical limit and compare with the spectrum 
( 4.3|) (n = N) deduced from the F-flatness conditions. In the next section, we shall 
use the matrix model to solve for the superpotential in each of the massive vacua 
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and this uncovers the action of the whole of the modular group. 



The eigenvalues of $ in the confining vacuum are given in (|3.28|) . (Recall that 



the eigenvalues in this vacuum have the property that they vanish in the classical 
limit.) The first thing to notice is that the modular parameter of the elliptic function 
parametrizing the eigenvalue spectrum is tr/N with tr = r — (iN/n) In A, where r 
is the gauge coupling. Hence S-duality, if realized, must act on the renormalized 
coupling tr as tr — > — 1/tr rather than r itself. Under this action, the modular 
parameter of the elliptic functions changes as tr/N — > — 1/{trN). But now we want 
to take the classical limit tr — > ioo and this can be done easily after performing a 
modular transformation —1/(trN) — > trN on the (quasi-)elliptic functions. Using 
the modular transformation properties of the theta functions (see Appendix A) we 
find 

fa) — Sx NM ^ e ~ tNTR " x + ^r^ ; x G ^ ] ■ (5 - 3) 

2A sm 2A sm | 

Apart form the additive constant, this coincides precisely with the expression Eq. ( |4.3[ ) 
in = N) for the eigenvalues in the Higgs vacuum of the deformed SU(iV) Leigh- 
Strassler theory with (3 replaced by t r [3. Equivalently, the Higgs and confining 
vacua of the theory with deformation parameter (3 are simply exchanged under the 
combined operation tr — > —1/tr and (3 — > (3/tr 

The overall factor of tr in the above expression is expected because the eigen- 
values have modular weight one. The fact that the additive constants are different 
is to be expected: even in the context of M = 1* theory the matrix model approach 
only yields results that are S'-duality covariant up to a vacuum independent shift [7]. 
This additive shift will be discussed in more detail in the next section. We should 
also point out that after S'-duality on the spectrum of eigenvalues in the confining 
vacuum we have taken a classical limit above. In this limit, the exact transformation 
properties of the couplings A, (3 are not visible. These will become clear in the next 
section. 



The results of this simple analysis are that there is indeed a form of S-duality 
that relates the Higgs and confining vacua of the theory with deformation parameter 
(3. In particular this duality acts on the bare couplings of the Leigh-Strassler theory 

as 

Tr — -1/TR , — (3/TR . (5.4) 

This is consistent with j3 having modular weight -1 as expected. Recall that the 
bare theory has 3 couplings r, A and (3 with one non-trivial relation, so that ( |5. | ) 
determines the action of S-duality on the bare theory up to the unknown relation. 
In the next section we expand the discussion of duality to include the whole of the 
modular group SL(2, Z) and deduce the action on all three of the bare couplings. 
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5.2 Multi-cut solutions of the matrix model and SL(2, Z) 



Now that we have described the structure of massive vacua of the Leigh-Strassler 
deformed theory, we can return to the matrix model and ask whether we can solve 
for the exact superpotential and condensates in each of these vacua. This will turn 
out to be very important because it will prove that there is an action of electric- 
magnetic duality on the Leigh-Strassler deformation of the M = 4 theory. In fact 
this extends to the action of the full SL(2, Z) modular group. 

The way to use the matrix model to solve for the theory in an arbitrary massive 
vacuum is described for the M = 1* theory in [8]. One finds that the same techniques 
can be used in the deformed matrix model described in Section 3. In order to describe 
the massive vacua with unbroken SU (N/p) symmetry, one needs a multi-cut solution 
of the matrix model. Of course the matrix model computes the superpotential in the 
U(iV) theory and we then use fl2.9p to deduce it in the SU(iV) theory. 



We start with the observation that the classical eigenvalues in the massive vacua 
with unbroken SU (N/p) are given by (|4.3| ) with n = p, each with degeneracy N/p. 
The only thing that will be important for us will be the range of the phases of 
these classical eigenvalues. The classical eigenvalues of the U(iV) gauge theory also 
have exactly the same phases in the corresponding massive vacua. As we have seen 
explicitly in the case of the confining vacuum, the eigenvalues of the field theory are 
essentially proportional to exp[z]. This directs us to consider a p-cut solution of the 
matrix model with the cuts in the z-plane situated at 

[a + i(3(p + 1 - 2j)/2, b + i(3(p + 1 - 2j)/2] ; j = 1, . . . , p (5.5) 

where these positions are dictated precisely by the phases of the classical eigenvalues. 
The classical eigenvalues are clustered in groups of N/p at each of these points, but 
the pairwise effective interactions induced in the matrix model will cause each cluster 
to spread out into a cut. 

Following [8], we then make an ansatz where the cuts have the same length and 
the same filling fractions. It can then be shown that G(z) defined in ( |3.8|) again 
only has two cuts, as in the confining vacuum, but where the cuts are shifted to 
±[a + iPp/2, b + ifip/2}. This modifies the map z(u) in (|3.13|) by a shift (3 — > p(3\ 

_ _J_ 01(01^(^2*1 -P/V4|f) 
P [ ) 2tanf 0UP/V2|f)0i(W2£i+P/V4|f) ' 1 J 



Similarly 



G «"» = CxlZ'i ■ w + aM ~ ^ lPPM] ■ (5 - 7) 
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When the Dijkgraaf-Vafa recipe is followed one finds that for each cut there is 
an associated filling fraction of eigenvalues Nj/N which translates to an effective 't 
Hooft coupling Sj = g s Nj with Yl Sj = S. For the massive vacua we note that the 
filling fractions in each of the p cuts are the same and equal to N/p = q 



2 ™ S i = ^jjr> 3 = h---,P, (5.* 



where h(r) is the same function that we encountered in ( [3.20] ) with (3 — > p/3: 



/zM* cos g W /2|f) 
[ ) 4A 2 sin 3 f e[(p(3/2\f) ■ ^ 

In addition one also needs the variation in the genus zero free energy upon trans- 
porting an eigenvalue from any given cut to infinity. After minor modifications this 
is also given by the same expressions that were encountered in the one-cut solution 
for the confining vacuum: 

±^-- f ^-m-^±s,+^ (5,0, 

Then the superpotential for the U(iV) theory at any of its massive vacua is 

wT } = ^±^-2*ir± Sj . (5.11) 

3=1 3 3=1 



Extremizing with respect to r gives 



(5-12) 



and therefore we find that the value of the quantum effective superpotential for the 
U(iV) theory is 



rV{N) _ N^M 2 cos f e l (pP/2\p 2 r R /N) iV/iM 2 
4pA 2 sin 3 f ' 9[(p(3/2\p 2 T R /N) 4A 2 sin 2 f 



In all of the above we should remember that p is a divisor of N. 



The result for SU(iV) follows by the application of Eq. 



w sv(N) _ pN^iM 2 O^pPWtr/N) NfiM 2 
eff ~2A 2 sin/? e 1 {pf3/2\p 2 r R /N) 4A 2 sin 2 ^' ' 



This gives the result in only one of q — N/p vacua of this type. The results in the 
remaining vacua follow by the replacement 

tr — >r R + k/p; fc = 0, -1 . (5.15) 
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The Higgs vacuum corresponds to p = N whilst the confining vacuum considered 
previously has p = 1 and k = 0. 



One interesting and immediate check on the form of the superpotential follows 
from the fact that it is manifestly divergent when the deformation parameter takes 
the special values 

Owk 

P= , keZ . (5.16) 

p 



But this matches precisely the values of (3 in Eq. (|4.4j) for which the classical eigen- 
values of $ for the massive vacua diverge (since they consist of q blocks of size n = p). 
Another check follows by taking the classical limit of the superpotential (|5.14 ). Using 
( |A.3c|) , one finds 



W3™ -~ <Tr - (V^ " l) . ("7) 

4X Z sin | V tan ^ / 

which matches the values of /i(Tr <f> 2 ) obtained in ( |4.5| ) from the F-flatness condition 
precisely. 

Now that we have the superpotential in each of the massive vacua, we can now 
investigate how the full SL(2, Z) duality group acts on the theory. In the Af = 
1* theory described in [7,8], we established that the matrix model yields a result 
for the superpotential which is only duality covariant when shifted by a vacuum 
independent constant. This kind of additive ambiguity is ubiquitous in this subject: 
see [7,9]. Naturally, the same kind of ambiguity will be present in the Leigh- Strassler 
deformation where we claim that the appropriate definition of a duality covariant 
superpotential is 

W s JW = W^ + C(t r ,(3,X), (5.18) 



where the vacuum-independent constant is given by either of the expressions 
° {TR ' /3 ' A) - iAWf " 2^m7 ' WWn) 12Psin/3 ^™ 



4A 2 sin 2 | vrA 2 sin/3 V n 



71 J 



(5.19) 



Here, the Weierstrass (^-function is defined relative to a torus with periods u)\ and 
0J2 and r R = U2/UJ1 and E 2 {t r ) is the 2 nd Eisenstein series (see the Appendix of [7] 
for definitions and references). One property that the C has is that it reduces to the 
correct form needed in the M = 1* theory in the limit /3 — > 0. 
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Our claim is that the duality-covariant definition of the superpotential in the 
each of the massive vacua labelled by (p, k), p ■ q = N and k — 0, . . . , q, is then 

W SU(N) _ NyM* / fl(gg/2jf) e[((3/2\r R ) ft x 
where 

f = (pr R + k)/q . (5.21) 



Using the modular properties of the theta functions and the 2 nd Eisenstein series, 
( |A.2a|) and (|A.5|) , one can easily show, for example, that 



~ sin 3 ~ 

W eS (ft/r R , -1/t r , A) = r R . -f ■ W eff (/3, r fl , A) . (5.22) 

Conf Sm{P/TR) Higgs 

Notice that the coupling /? transforms with unit weight under the modular group. 
Note also the rather unconventional transformation of the superpotential itself: this 
we will have to cure by an appropriate definition of the transformation of the coupling 
A. 

A consideration of the other massive vacua show that they all lie on a single 
orbit of the SL(2,Z) duality group. For the element 

„= ^ ^ e SL(2,Z) , (5.23) 

which relates two vacua A = cr{B), we have 9 

_§_ OTR+b 



sin P 



sm[f3/(cr R + d)] 



B 



From this we can deduce the action of o G SL(2, Z) on the bare couplings of the 
theory; clearly 

ar R + b p 

°(t r ) = — , <r{p) = — , (5.25) 

ct r + a ct r + a 

confirming that f3 transforms with modular weight -1. However, we can say more: in 

order that the superpotential transforms with a definite modular weight, so that it is 

consistent with the group action, we can deduce the transformation of the coupling 

A. Recall, that on the marginal surface of the Leigh-Strassler deformation A is a 

function — albeit unknown — of the other two couplings and so must transform. In 

addition, since in the limit /3 — > 0, the superpotential transforms with modular 

weight 2, we expect this to be preserved for /3 ^ 0. This requires the combination 

A 2 sin ft that appears as an overall factor in the superpotential must transform with 

unit modular weight: 

(X (A 2 sin/3) = ^^. (5.26) 
v ' cr R + d K ' 



3 Thc actual action of a on a vacuum with labels (p, k) is not difficult to elucidate. 
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6. An exact elliptic superpotential 

Another method that has been used to calculate the exact values of the superpotential 
in all the vacua of the M = 1* theory is to compactify it on a circle of finite radius [9]. 
The effective superpotential is then a function of the dual photons and Wilson lines 
of the abelian subgroup U^l)^ -1 C SU(iV). These comprise N — 1 complex scalar 
fields X a , a = . . . , N (with ^2 a=t Xa = 0) which naturally live on a torus of complex 
structure r because of the periodicity of each dual photon and Wilson line. The 
superpotential describing the M = 1* deformation is therefore constrained to be an 
elliptic function of the complex scalars X a . In [9] the exact answer was found to be 

where p(z) is the Weierstrass function defined on a torus with half periods uj\^ 
with t = U2/0J1. What is particularly useful about this superpotential is that it 
is independent of the compactification radius, and therefore, yields results that are 
valid in the four-dimensional limit, and it encodes all the vacua of the M = 1* theory: 
both massive and massless. 

The question is whether there is an elliptic superpotential that describes the 
Leigh-Strassler deformation? Similar arguments to those used in [9] suggest that 
there should be, but are not powerful enough to determine the exact answer. How- 
ever, there is an obvious candidate for appropriate marginal deformation of ( |6.1| ), 
which turns out to be correct. Specifically, we can replace the Weierstrass function 
p(z), by its unique deformation in the space of even elliptic functions. 10 In particular, 
the Weierstrass function has one double pole at the origin and we can consider the 
elliptic function of order two obtained by splitting the double pole into two single 
poles. Remarkably, a superpotential of precisely this form reproduces our results 
from the matrix model for the superpotential in each of the massive vacua. This 
superpotential can be written as, 

f a^b 

The Weierstrass-zeta function ((u) can be defined via p(u) = —('( u ) where <p{u) is 
the Weierstrass function for the torus with complex structure tr. Then in the limit 
(3 = this reproduces (|6.1| ). What is remarkable is that the massive vacua of the 
M = 1* theory for which the X a lie on sub-lattices are also critical points of (|6.2j) 
since this only relies on the fact that the two-body interaction is an even elliptic 

10 The function must be even to respect the Weyl group of SU(7V) 
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function. For instance, for the p ■ q = N massive vacuum with fc = 0we have 

v q 

Furthermore, the values of the superpotential are precisely given by ( J5.20 ). 



X a e{^ + ^, r = 0,...,p-l, S = 0,...,g-l). (6.3) 



A very strong check on the form of the superpotential we have conjectured in 
( |6.2|) is possible because of the existence of new vacua which are not present in 
the limit (3 = 0. For example, it is easy to show that X a = 0, Va, is a critical 
point. Actually this is an SL(2, Z)-invariant vacuum. Now we can take the classical 



limit, t — > zoo, after having shifted by minus the vacuum-independent shift ( |5.19| ), 
and compare with the analysis of the vacuum structure in Section 4. In the case 
of SU(3), the classical limit of the superpotential W^xs 1 — C in the vacuum with 
X a = precisely equals the value calculated from the new vacuum solution of the 
F-flatness conditions (f4.9|). We expect this to generalize to SU(A^), however, the 
analysis remains to be done. 



The properties and construction of this elliptic superpotential will be described 
in a separate publication; however, we mention the fact that for M = 1*, the el- 
liptic superpotential is, more or less, the (complexified) Hamiltonian of the elliptic 
Calogero-Moser integrable system which described N particles interacting though 
pair-wise forces proportion to p(X a — X b ). The intriguing question is whether this 
system can be deformed 

p(X a - X b ) -> ((X a -X b + M/tt) - ((X a -X b - Pux/n) (6.4) 

whilst preserving integrability. 

Finally we comment on a remarkable correspondence between our results for 
the Leigh- Strassler deformation and the exact superpotential of a certain relevant 
deformation of the Ai-quiver theory studied in [14]. Specifically this is the theory of 
an SU(iV) x SU(iV) gauge multiplet with gauge couplings g\ and g 2 , with two chiral 
multiplets A and B in the (N, N) representation and two chiral multiplets A and 
B in the (N,N). The theory also contains two chiral multiplets $1 and $2 in the 
(adj,0) and (0, adj) representations respectively with classical superpotential 

W = Tr{{g 1 Qi + M/2)(AA + BB) + (g 2 $2 + M/2)(AA + BB)} . (6.5) 

This theory actually has M = 2 supersymmetry but we further deform it by adding 
the operator, 

/zTr _ $2) . ( 6 . 6 ) 

to the classical superpotential which preserves M = 1 supersymmetry. This theory 
is a relevant deformation of the M = 2 superconformal theory obtained by setting 
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\i — M — 0. The exact vacuum structure and superpotential was determined in [14]. 
The theory has N J2 d \ N d massive vacua. Remarkably, in J2 d \ N d of these vacua, we 
find exactly the same value of the superpotential as in each of the massive vacua of the 
SU(iV) Leigh- Strassler theory studied in this paper. Specifically, this agreement holds 
if we identify the parameter (3 of the Leigh- Strassler deformation and the complexified 
coupling of one of the SU(iV) factors in the quiver theory n = Ajri/gf + 9\/2ir (or 
z = 2z7TTi to use the notation of [14]) : 

, = !™ . (6.7) 

7T 

In the brane set-up for the quiver model z is the relative separation of the two NS5- 
branes. We also identify the parameters M and \x of the quiver superpotential with 
parameters of the same name in the deformed Leigh-Strassler superpotential ( |1 . 1| ) . 

This connection suggests that there should be an RG flow from the superconfor- 
mal Ai-quiver to the Leigh-Strassler fixed line. At least in one special case, such a 
flow is known to exist. In [23], Klebanov and Witten considered the quiver theory 
in the case of equal gauge couplings g\ = #2 and zero masses for the bi-fundamental 
multiplets, M = 0. In this case one can integrate out the adjoint chiral fields $i 
and $2 and induce a marginal quartic superpotential for the bi-fundamentals which 
corresponds to D3-branes placed at a resolved conifold singularity. Moving onto the 
Higgs branch where SU(iV) x SU(iV) is broken to the diagonal SU(iV), we find an 
effective theory of three chiral superfields transforming in the adjoint of the unbroken 
SU(iV) which coincides with M = 4 SUSY Yang-Mills. Repeating this exercise with 
unequal gauge couplings g\ 7^ we find that the Leigh-Strassler marginal defor- 
mation as well as other irrelevant operators are induced. However, it is not obvious 
whether this flow is related to the exact agreement we have found above. 
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Appendix A: Some properties of elliptic functions 

The (quasi-) elliptic functions that we need, p(z), ((z) and 9i(z\r), are all stan- 
dard and their definitions and properties can be found in [20]. In our conventions, 
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the basic torus has periods 2oj\ and 2oj 2 with complex structure r = lj 2 /ui. We also 
define q = exp(?7rr). Below we list some properties of these functions that will be 
useful. 



An important relation between the Weierstrass and Jacobi theta functions is 

C(M) _ m u = jL«jMpm . (a.i) 

The Jacobi theta functions have the following modular properties that will be useful 

9 1 (x\t) = z(-zr)~2 exp(-z:r 2 /W)#i(£/V| - 1/r) , (A.2a) 
9 3 (x\t) = (-ir)~5 exp(-ix 2 /7rr)0 3 (a;/r| - 1/r) . (A.2b) 

In order to take the classical limit, we will need the expansion 

oo 

0i(z|t) = 2 J2{-l) n q [n+1/2)2 sin(2n + l)x , (A.3a) 

n=0 

oo 

3 (<e|t) = 1 + 2^(-l)V 2 cos2nx , (A.3b) 

n=0 

#'(x|t) °° a 2n 

lV ( =cotx + 4y^ - - sin2nx . (A.3c) 
9\{x r) 1 — (r n 

v 1 y n=l 

We also need the 2 nd Eisenstein series [22] 

E 2 (t) = 1 -24 J2°i(n)q 2n , (A.4) 

n=l 

where 0i(n) is a sum over each positive integral divisor of n. It has the following 
non-trivial modular transformation 

fir 

E 2 (-1/t)=t 2 E 2 (t) + — . (A.5) 
Appendix B: Calculating condensates 



In this appendix we show how the integrals of the type ( |3.25a| ) and ( |3.25b| ) can 



be computed to yield the condensates in the U(A) and SU(A) gauge theories. Let us 
first consider the function f(x) in terms of which the eigenvalues have been defined 

f(x) _ e,(nx-P/4\f) _ 9 3 (nx-n/2-P/4\f) _ 

J[ ) 9 4 (irx + (3/A\f) ~ 3 (7rx-7r/2 + /3/4|f)' L ' 1 ' 1 } 
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Clearly, the eigenvalues can be written either in terms of #3 or 64 since the only 
difference is in the range of the argument x. We also define the function g(x) 



9 1 (nx-i3/i\f) 



0i(7rz + /3/4|f) ' 
The functions / and g have some nice properties 



(B.2) 



g(x + l)=g{x); g{x + f/2) = e^ 2 f(x) ; g(x - f/2) = e^ 2 f(x) , (B.3) 

which follow from standard identities. Note also that the function g(x) has a simple 
pole in its fundamental period at x — l — /3/4n. Now it is a simple matter to convince 
oneself that because of the periodicity properties of g(x) the contour integral around 
this pole and hence the residue there can also be written as 



Res g{x) 



Thus we see that 



X = -/3/47T 



f(x) dx 



-2i sin § 



f(x) dx 



1 giQg/2^) 
sinf 0[(O\f) 



(B.4) 



(B.5) 



From this key result it follows using Eq. ( |3.24| ) 



Tr$ 



N / dxcj)(x) 



NM 



2A sin 



I* 



1 W/gjf) 
taa|^09/2|f) 



1 



2Asin f W V{N) 



Mji 

(B.6) 

(Note that because <p(x) is a periodic function under x — > x + 1 the actual limits 
of integration can always be changed from (—1/2,1/2) to (0,1).) Similarly it also 
follows that in the SU(iV) theory 



Tr$) = iV / 0(x)dx = 



(B.7) 



To compute the quadratic condensates we need the integrals f 2 (x)dx. These 
can also be calculated using the same trick as above and we find that 



f 2 {x) dx 



1 



2i sin (3 



Res g 2 (x) 



20' 1 (/3/2|f)0 1 (/?/2|f) 



-/3/4-7T 



sin^0i(O|f 



(B.8) 



Using the results Eqs. (B.5| ) and (|B.8|) we can readily compute the quadratic 
condensates. For the U(iV) theory we find 



(Tr$ 2 ) = N <f) 2 {x)dx 



NM 2 



NM 2 



2 A 2 sin/? 6[([3/2\f) 4A 2 sin 



2 p 



l -w. 



U(AT) 
eff 



(B.9) 
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For the SU(iV) theory we find similarly 




NM 2 #i 03/2 |f) NM 2 



-W^ (N) (B.10) 



2A 2 sin/5'^(/5/2|f) 4A 2 sin 2 f 



Similarly higher condensates could be computed explicitly if required and will in 
general depend on the residue of g n (x) at the pole x = —f3 /Arc 
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